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ABSTRACT 

We investigate the microlensing effects on a source star surrounded by a circumstellar disk, as a 
function of wavelength. The microlensing light curve of the system encodes the geometry and surface 
brightness profile of the disk. In the mid- and far-infrared, the emission of the system is dominated 
by the thermal emission from the cold dusty disk. For a system located at the Galactic center, we 
find typical magnifications to be of order 10-20% or higher, depending on the disk surface brightness 
profile, and the event lasts over one year. At around 20 ^va, where the emission for the star and the 
disk are comparable, the difference in the emission areas results in a chromatic microlensing event. 
Finally, in the near-infrared and visible, where the emission of the star dominates, the fraction of star 
light directly reflected by the disk slightly modifies the light curve of the system which is no longer 
that of a point source. In each case, the corresponding light curve can be used to probe some of the 
disk properties. A fraction of 1G~^ to 10~^ optical microlensing events are expected to be associated 
with circumstellar disk systems. We show that the lensing signal of the disk can be detected with 
sparse follow-up observations of the next generation space telescopes. While direct imaging studies of 
circumstellar disks are limited to the solar neighborhood, this microlensing technique can probe very 
distant disk systems living in various environments and has the potential to reveal a larger diversity 
of circumstellar disks. 

Subject headings: gravitational lensing — circumstellar matter 



1. INTRODUCTION 

During the last decade, gravitational microlens- 
ing has become a powerful tool for analyzing 
the distribution and nature of compact dark ob- 
jects in the halos of our and nearby g a laxies . 
The idea was first suggested by IPaczvhskil l)198fif) 
and nowadays a number of monitoring programs 
are ongoing and have detected several thousands 
of m i crolensing event candidates l|Al999 k ei^ ^1.1 
Tflflfil lAnbonrg et a,lJ ITOfl^ ITTdalski et alJ II 99,1 
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20021 iBozza et al.l 1200^ . In addition, several authors 
have pointed out that microlensing events themselves 
can yield useful information about the sources, i.e., 
the lensed st ars, by taking into account their angular 
extent (e.g., iGouldl 1^9 4: Simmons. Ncwsam, & Willis 
Valls-Gabau d ,.1998.: .Hendry et al.. 
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In this paper, we discuss another microlensing appli- 
cation: the magnification of a star surrounded by a cir- 
cumstellar disk. The microlensing properties of such a 
system appear to be quite interesting: in the optical 
and near infrared bands, the emission is dominated by 
the star light and classical .Paczvh ski light curves can be 
observed, whereas in the mid/far-infrared, the emission 
from the cold dust of the disk dominates and produces 
more complex light curves. Microlensing translates spa- 
tial structures into temporal ones and allows us to "re- 
solve" the disk. Combining information from both wave- 
length ranges can place constraints on some of the disk 
properties. In addition, because of the presence of a disk, 
a fraction of the star light is refiected by the dusty parti- 
cles. As this light contribution no longer originates from 
a point-like source, it will cause the microlensing light 
curve to slightly deviate from that for star-only in the 
optical and near-infrared. This effect can also be used to 
constrain disk properties. Direct imaging studies of disks 
are possible but are limited, in the foreseeable future, to 
the vicinity of the solar system. The microlensing tech- 
nique allows us to access systems that are much further 
away, i.e., up to ^ 10 kpc. This technique can probe cir- 
cumstellar disks in various environments with different 
metallicities and has the potential to reveal the diversity 
of such systems. 

The characteristics of circumstellar disk microlensing 
light curves differ from those of stars: the total magnifi- 
cation of the disk in the mid-infrared is expected to be 
significantly smaller than that of a star in the optical. In 
addition the relative duration of the microlensing event 
for the disk is expected to be much longer. Indeed, the 
typical radius of a circumstellar disk a is of the order of 
a few tens of AU, t hough it depends on the age and type 
of the central star l|Zuckermanl2fl0ll) . For a source at the 
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Galactic center and a lens of 1 Mq located halfway, the 
Einstein ring radius te in the source plane is about 8AU, 
which is much smaller than the source (disk) size. We 
therefore have an extreme finite source effect and only 
that part of the source inside the Einstein ring is effi- 
ciently magnified. For a large uniform face-on disk, the 
expected overall fracti onal flux exc ess caused by the lens- 
ing effect is 2{f-E,/af l)Agolll200^ . which is at the level 
of ~ 10% if the disk radius is a few times larger than 
te- The microlensing timescales are also different. For a 
star, the duration of the lensing event is, i-e., the time 
required to cross the Einstein radius, is about one month. 
For a circumstellar disk, this timescale is about {a/rE)tE, 
which ranges between half a year to a couple of years. If 
a microlensing event of a (young) star is discovered (in 
optical or infrared), one can then (sparsely) monitor the 
object in the mid-infrared and look for the excess charac- 
teristic of circumstellar disks. If a disk can be detected, 
then observing its flux as a function of time will pro- 
vide us with constraints on the disk size and brightness 
distribution as the lens differentially resolves the disk. 

In this paper we show that next generation space tele- 
scopes will be able to detect microlensing events of cir- 
cumstellar disks at the Galactic center. The outline of 
this paper is as follows. In § 2, we present some charac- 
teristic parameters of circumstellar disks. In § 3, we in- 
troduce the formalism describing the observed light curve 
of a disk. In § 4, we show how the disk parameters are 
related to the shape of the light curve and the behaviors 
of light curves in different observational bands. We then 
investigate the detectability and the occurrence of this 
phenomenon in § 5. Finally, we summarize and discuss 
our results in § 6. 

2. CHARACTERISTICS OF CIRCUMSTELLAR DISKS 

Circumstellar disks evolve through two main stages: 
during the first phase, the protoplanetary disk composed 
of gas and dust is optically thick and provides a parti- 
cle reservoir for planet formation. By the time stars of 
~ 1 M0 reach ages of ^ 1 Myr, the accretion disk sur- 
rounding the stars has lost its gas, becomes optically 
thin and readily accessible to observation at many wave- 
lengths. It is referred to as a debris disk. The dust 
particles in the disk absorb a fraction of the star's pho- 
tons and re-emit them at longer wavelengths. The pa- 
rameters of some debris disk s observed in the vic inity of 
the solar system are listed in lZuckerman fc Sond l|200H) . 
Current observations show that typical systems have a 
star temperature of about 6000 K and the excess infrared 
emission from the disk can be well fitted by a blackbody 
spectrum with an effective temperature of about 80 K. 
This does not imply that the disk is isothermal but is 
simply a useful way for describing the disk spectral en- 
ergy distribution. 

Apart from the re-emitted light, some of the star pho- 
tons intercepted by the disk are reflected. The relative 
amount of light being reflected and that being absorbed 
and re-emitted depends on the albedo of the disk, which 
is of the order of a few tens of percent. For simplicity, we 
assume the albedo to be 0.5 so that the disk seen through 
the reflected light and through the re-emitted light has 
the same luminosity, which we denote as a fraction r of 
the star's bolometric luminosity. 

It has been empirically found that the fraction r de- 



pends on the age of the syst em, r (x age~^'^^ for ages 
between 10*'-^ and 10^-^ years l|Zuckerman fc SontJ2 004). 
For a 10 Myr system and a 100 Myr system, r ~ 10~^ 
and ^ 2 X lO""*^, respectively. The expected spectral en- 
ergy distribution of such systems located at the Galactic 
center (a distance of 13 = 8 kpc), is illustrated in Fig- 
ure ^ for different values of t. In this plot we have not 
considered extinction effects as they are in general small 
for the wavelength ranges we consider below. As can be 
seen, the infrared excess caused by the thermal emission 
of the disk decreases rapidly as the system ages. In the 
figure we have assumed that the reflected light from the 
disk has the same spectral energy distribution as the star. 

Circumstellar disks are expected to have an inner 
disk radius set by the su blimation temperature of 
the dust grain l)Natta et al. 2001) with radii rang- 
ing from around one to a few AU, and in addition, 
theoretical models predict that planets in disks can 
inter act with the disk material and open wide radial 
gajjs ifPapaloizou & LinI 179841 iPaardekooner fc Mellemal 
l2004{l . once they achieve sufficient mass. The width of 
the gap is expected to be ~20% of the orbital radius, 
which would decrease the disk surface density in a 1 
AU annulus at 5 AU for example. In addition disk 
structures such as rings, warps and blobs are expected 
as a result of the presence of planets llMoran et alJ 
2004 'Kuchncr & Holman 2003; Wilncr ct al. 200^ 
Wvatt fc Dent, .2002; .Quillen fc Thorndike ,200^ 
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imprinted in the microlensing light curve of the system 
and, in some cases, might enhance the detectability of 
the disk. 

At the Galactic center, the angular size of a 30 AU 
disk is only ~ 0.004". Spatially resolving such system 
in the visible and mid-infrared would require a telescope 
larger than 30 m, and 1 km respectively. With microlens- 
ing, however, the spatial information is translated into a 
temporal one. As shown below, measuring the microlens- 
ing light curve allows us to probe the disk geometry and 
surface brightness profile. 

3. THE FORMALISM OF MICROLENSING OF DISKS 

The characteristic scale for microlensing is the Ein- 
stein radius of the lens. For our purpose in this paper, 
the Einstein radius is defined in the source plane, which 
reads 

1/2 
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where M is the lens mass, x — Di/Dg, and Dg and Di 
are the distances to the source and lens, respectively. For 
a source at the Galactic center and a 1 Mq lens located 
half way, x — 1/2 and te — 8 AU. The corresponding 
time scale for the source to cross the Einstein radius is 
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Fig. 1. — Spectral energy distribution of a star with T = 5800/f 
surrounded by a dusty debris disk with T = 80 K. This system 
is put at the Galactic center (distance D = 8 kpc). The reflected 
(dashed curves) and re-emitted (dotted curves) lights from the disk 
are calculated for three values of the disk-to-star luminosity ratio: 
T = 10-2, 10-3, and 10-*. 



where v = DsfJ-rci and /irci is the relative proper motion 
between the source and the lens. 

Considering the central star as a point source, the mag- 
nification A of the source can be shown to depend only 
on the projected separation u of the lens and source, in 
units of te of the lens (^Pacz vuski 198^, namely 



Aiu) 



u \Ju?- -I- 4 



(3) 



The relative motion between the source and the lens 
changes u and therefore the magnification varies as a 
function of time. Microlensing events of point sources 
can be identified from this characteristic light curve. 

^Hcvrovsky & Locb (1997) considered the microlensing 
of inclined accretion disks of active galactic nuclei. In 
this case, the Einstein ring radius is much larger than 
the size of the source and large magnification values can 
be reached. The case considered in this paper is dif- 
ferent: the size of Galactic circumstellar disks is larger 
than the Einstein ring radius leading to lower magnifica- 
tions and longer timescales. For the magnification of very 
large sources, a number of analyti cal approximations 
have been derived in the literature ffle vrovskv fc Loebl 
[1997, AeoI 2003). In our case, the Einstein ring is not 
necessarily much smaller than the spatial extents of fea- 
tures in the disk, therefore we proceed with the exact 
calculations. 

In this paper we do not intend to model the disk in 
great details. For estimating the lensing signal we will 
use different surface brightness profiles in order to repre- 
sent the changes in dust and temperature distributions. 
Let us now consider an axially symmetric disk with sur- 
face brightness profile /(a), where a is the radial distance 
to the central star. The projection of the disk onto the 
sky plane is an ellipse. The configuration of the disk is 
fixed given its center position (i?, $) with respect to the 
lens, the inclination angle i, and the orientation angle 




Fig. 2. — Geometry of lensing of a circumstellar disk. The dotted 
ring around the lens indicates the Einstein ring. The circumstellar 
disk is assume to have an hole inside an inner radius and a cutoff 
outside an outer radius. Projected onto the sky plane, the disk 
appears to be an ellipse with the orientation angle of the major 
axis being cAq. Any point on the disk is specifled by its coordinates 
(r, 0) in a polar system centered at the star. The angles and 
4i are with respect to the direction of the relative motion v. The 
position of the star (or the disk center) is specified by (i?, #) with 
respect to the lens (<I> is relative to the direction of v). The impact 
parameter of the star is denoted as . 



of the projected major axis (with respect to the direction 
of the relative motion). The geometry and meanings of 
symbols are shown in Figured In all our discussions be- 
low, quantities with length dimension are in units of the 
Einstein ring radius te and surface brightnesses are in 
units of flux/f|. We consider a geometrically thin disk 
with isotropic emission at every point. With no lensing 
effect, the total flux of the disk can be calculated without 
projection and is given by 



Fn = 



f{a)adad(j) — 2tt f{a)ada 



(4) 



If the disk is lensed, the total fiux is the unlensed fiux 
convolved with the magnification profile A{u), 



F = 



A{u)fp{r, (pydrdcj), 



(5) 



where fp is the projected surface brightness profile, u = 
[R"^ + 2Rrcos{^ — 0) + r^]^/^ is the distance from the 
lens to a (projected) disk element. Here (r, 0) are po- 
lar coordinates with origin at the center of the star (see 
Figure [21 • The integration becomes easier if we change 
the variable r to the radius a along the major axis ac- 
cording to the relation a — 



'1 + sin^((?!)- 0o)tan2i 

The projected surface brightness fp is simply f{a)/cosi 
Equation ((Sj then has the form 



F = 



A{u) 



cosi 1 + sin (0 — 0o) tan i 



dadcj). (6) 



The overall magnification of the disk fiux is F/Fq. The 
temporal dependence of the magnification is through the 
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position of the disk center {R, $) with respect to the lens, 
R{t) = J t2 + ul and = arccos^^ ^ (7) 

where uq is the impact parameter in unit of te, and t 
is the time in unit of is = ^e/w defined in equation ^ . 
We set i = at $ = tt/2. 

From equation 0, we see that the hght curve of the 
disk depends on the magnification kernel A{u) and the 
geometry and the surface brightness profile of the disk. It 
should be noted that A{u) can be directly inferred from 
observations of the light curve of the star in optical or 
near-infrared band. The disk light curve can be used to 
constrain /(a) and the geometric parameters only. 

Before investigating the expected light curves de- 
scribed by the exact results of equation lO, it is in- 
teresting to get a more intuitive understanding of this 
equation by considering the following approximation: if 
the spatial features of the disk are much larger than the 
Einstein ring radius, the calculation of the magnification 
can be simplified by approximating A{u) — 1 as 27rj£)(u), 
where Sd is the Dirac-(5 function. The increase in flux 
AF = F — Fq caused by the lensing simply traces the 
local surface brightness of the disk at the position of the 
lens (i.e., r = — R). The result is 

AF = 27r/p(r,(/)) = 27r/(a)/cosi, (8) 

where r = R{t), (p = ^{t) — tt, and a is related to r and 
(j) as mentioned above equation ©. The excess flux is 
just twice the Einstein ring area times the surface bright- 
ness at the lens pos i tion, w hich is the same as derived by 
IHevrovskv fc Loebl JTool ( s ee the ir eq.[34]) for an el- 
liptical source and by lAgol ll2?ffll (see his eq.[10]) for 
general large sources. In fact, the result of_the mag- 
nification for a large source presented in lAgoll l)2003f) 
(his eq.[4]) can be derived as above by noticing that 
A{u) — 1 ~ 2tt6d{u). Approximate calculations of the 
magnification at the edge of a large source can be found 
in Agol ( 2003). 

It should be noted that the light curve alone is not 
sufficient to determine the physical size of a disk and 
what we infer are all in units of the Einstein ring radius. 
Additional information on the size of the Einstein ring 
radius could be obtained from the star-only light curve 
if the finite source effect and the parallax effect are de- 
tected. The former gives the angular size of the Einstein 
ring radius with respect to that of the source star (e.g., 
IGouldll994j) and the latter measures the Einstein ring ra- 
dius projected ont o the observer's plane ( see, e.g.. lGouIdl 
119921 iJiang eral]l200l ICould et al.l2"oo|l . Since the im- 
pact parameter needs to be small for a measurable finite 
source effect, it is unlikely that the lensing of a system 
with circumstellar disk can measure the angular size of 
the Einstein ring. Even if we have both the relative an- 
gular size of the Einstein ring and its projected size onto 
the observer's plane, we still need the source distance 
in order to infer the size of the Einstein ring projected 
onto the source plane. The radial velocity of the source 
star could provide some information on the source dis- 
tance. If the blended light in the microlensing event is 
caused by the lens, which can be verified astrometrically, 
information about the m ass and distance of the lens can 

be obtained l)Ghosh et al 2004) . Finally, if we observe 

the event at different wavelengths, the temperature as 



a function of radius in units of the Einstein ring radius 
can be constrained. Then, by knowing the dust particle 
albedo, measuring the central star temperature, and as- 
suming thermal equilibrium, the corresponding physical 
radii can be inferred. In general, we may have to rely 
on some statistical priors to constrain the physical size 
of the disk. 

4. OBSERVATIONAL SIGNATURES 
4.1. Mid- and Far-infrared: Re-emitted Light 

In the mid- and far-infrared, the luminosity of the sys- 
tem is dominated by the light re-emitted by the disk and 
the contribution of the star can be neglected. The light 
curve of the disk depends on its geometrical properties 
(the size, the inclination angle i, the orientation angle 
00 of the projected major axis, and the impact param- 
eter Uq) and its surface brightness profile. To show its 
dependence on various parameters, we assume that the 
disk has an inner radius and an outer radius ao with a 
power-law surface brightness profile /(a) oc a" between 
these two radii. The light curves are shown in Figure |31 
with Ui, ao, i, 4>o, uq, and a varied, one at a time. 

As shown in panel (a), the existence of an inner ra- 
dius of the disk causes a corresponding dent in the light 
curve because the disk is less efficiently magnified as the 
lens comes close to the inner radius. The dent becomes 
deeper as the inner radius becomes larger. The width of 
the dent is proportional to the inner radius. Similarly, 
radial gaps opened by planets in the disk will also exhibit 
such dents in the light curve. Panel (b) shows cases for 
which the outer radius varies: the smaller the disk, the 
larger the magnification as the percentage that can be 
effectively lensed becomes larger. The overall magnifica- 
tion is inversely proportional to — af ^ a^, and the 
duration time is proportional to ao- 

The area of the disk projected onto the sky plane de- 
pends on the inclination angle, i.e., proportional to cosi. 
With other parameters fixed, a more inclined disk would 
have a higher overall magnification, as shown in panel 
(c), because its projected area is smaller and it becomes 
more efficiently magnified. On the other hand, the ori- 
entation of the projected disk does not change the area, 
and has therefore little effect on the maximum magnifica- 
tion. This is shown in panel (d). However, the duration 
has a dependence on the orientation, since it changes 
the path length that the lens crosses the disk. Except 
for cases with (j)o = 0° and 90°, the orientation of the 
disk (plus the non-zero impact parameter) also causes 
the light curve to be asymmetric about t=0. 

In panel (e), we investigate the effect of the surface 
brightness profile of the disk. With a more concen- 
trated surface brightness profile, the disk appears to be 
less extended, and the maximum magnification increases. 
In our example, the maximum amplitude changes from 
^--^15% for a uniform disk to ~40% for a disk with surface 
brightness oc a~^. Slopes at the rising and fading parts 
of the light curve also probe the surface brightness pro- 
file, in the sense that a steeper light curve corresponds 
to a steeper brightness profile. As we mention in § |21 if 
the disk is much larger than the Einstein ring, the light 
curve simply traces the surface brightness profile along 
the lens trajectory. 

For a point source, the magnification is determined by 
the impact parameter uq- However, for a disk close to 
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Fig. 3. — Lensing light curves of the disk. In each panel, one 
disk parameter is varied as indicates at the upper-right corner with 
others fixed at values listed at the upper-left corner. The y-axis is 
the magnification. The disk parameters include the inner radius a; , 
the outer radius ao, the inclination angle i, the orientation angle 0o 
of the projected major axis, the surface brightness profile index a, 
and the impact parameter uq. Parameters with length dimension 
are in units of the Einstein ring radius. 



uniform, tlic impact parameter has little effect on the 
magnification of the disk flux, since the magnification is 
largely determined by the area ratio of the Einstein ring 
to the projected disk. Panel (f) shows that, because of 
the inner hole in the disk, a smaller impact parameter 
leads to a deeper dent in the light curve. 
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Fig. 4. — Lensing light curves with similar magnification. In the 
left panel, values of the outer radius ao and the inclination angle i 
are changed in a way to keep the projected area of the disk fixed 
[(a^ — ) cos i = constant]. In the right panel, additional changes 
of the orientation angle 0o are made to match the duration time. 

After having seen how the light curve responds to dif- 
ferent parameter changes, it should be noted that a num- 
ber of parameter combinations could be degenerate. To 
first order, the light curve of a uniform disk can be char- 
acterized by the maximum magnification and the dura- 
tion time. As we have discussed, the magnification is 
largely determined by the area ratio of the Einstein ring 
to the projected disk. For example, a larger disk that 
is highly inclined could have the same projected area 
as a smaller disk that is slightly inclined. The magni- 
fications for the two cases would be similar. The de- 
generate direction of the magnification roughly follows 
(Oq — of) cos i = constant. The left panel of Figure 01 



shows examples with (a^ — af)cosi fixed. As expected, 
the amplitudes of the light curves are similar. The du- 
ration time (in unit of t^) of the event is the length (in 
unit of te) of the chord the lens crosses. For a larger 
disk, it becomes longer, but this can change if we vary 
the orientation of the projected disk as seen in panel 
(d) of Figure 121 By adjusting the orientation parameter 
(jjo, we can have light curves that mimic each other, as 
shown in the right panel of Figure ^ The asymmetry 
of the light curves can break the degeneracy. But if the 
data are noisy, the solution would be partially degener- 
ate with respect to the size of the disk, the inclination 
angle, and the orientation angle. 

The above degeneracy can be broken by making use 
of the star-only light curve (in optical or near-infrared 
band). The orientation of the disk leads to a difference in 
the time of the center of the disk-only event with respect 
to the star-only event. This time difference, together 
with the impact parameter inferred from the star-only 
light curve, gives the orientation. Once the orientation 
is known, the duration time, which corresponds to the 
chord length, and the magnification, which is related to 
the area of the projected disk, can be used to infer the 
size and the inclination of the disk. For a realistic, non- 
uniform disk, the surface brightness profile has the effect 
to make the light curve more asymmetric if the orienta- 
tion is neither 0° nor 90°. This also helps to break the 
degeneracy between the orientation and the inclination 
(with the assumption of a circular, azimuthally symmet- 
ric disk in projection). 

4.2. The ^20fj,m Transition 

For a large range of disk-to-star luminosity ratios and 
temperatures, the flux of the star becomes as important 
as that of the disk in the ~ 10 — 20/xm wavelength range. 
As a point source, the star has its flux highly magnified, 
while the amplification of the disk flux is, as seen above, 
at the 10-20% level and lasts over a longer period. An 
interesting phenomenon in this case is that the lensing 
event appears to be chromatic if we compare the light 
curve in optical and that at 20/im. The reason is 
as follows: in optical, the baseline flux (unlensed flux) 
is from the star and the whole flux is magnified by the 
lens. At ^ 20/im, the baseline flux is the sum of those 
from the star and the disk, but only the star's flux is 
highly magnified by the lens. Therefore, the apparent 
magnification of the star-|-disk system is much less than 
that for the star. That is, if the fraction of the unlensed 
fluxes from the star and the disk are x and I — x, the 
magnification of the system is then xA^, + (1 — x)Ad, 
where the magnification A^, of the star is much larger 
than Ad, that of the disk. Figure El shows an example of 
the light curve for the system in optical and at ~ 20/im 
with the assumption that the star and the disk have the 
same flux {x = 0.5) at ~ 20/im. The apparent chromatic 
effect can happen if the system to be lensed has different 
spatial structures in two bands or there is a color gradient 
across the (large) system. 

4.3. Visible and Near-infrared: Reflected Light 

It is interesting to note that constraints on the disk 
can also be obtained at shorter wavelengths where the 
star flux largely dominates: a fraction t of the star light 
is directly reflected by the circumstellar dust (so strictly 
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Fig. 5. — Contributions of star and disk to the lensing signal. 
The plot assumes that the star and the disk have the same unlensed 
flux at ~ 20/xm. The dot-dashed curve shows the total Icnsed flux 
at ~ 20^im, which is the sum of that from the star (dashed curve) 
and that from the disk (dotted curve). The solid curve is the 
lensed flux in optical where only the star is visible. The lensing of 
the star+disk system is apparently chromatic because the system 
has different spatial structures in optical and in mid-infrared. 
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Fig. 6. — The top panel shows the classical microlensing light 
curve of a point source (dotted curve) and that of a star plus the 
disk in reflected light (solid curve). A fraction of r = 10~^ of the 
star light is assumed to be reflected by the disk. The bottom panel 
shows the relative flux difference: T{Aii/Ai, — 1). Parameters of 
the disk are the same as those used for the solid curve in panel (a) 
of Fig. El 



speaking, the star-only case we mentioned before is not 
exactly so). The reflected light has similar spectral en- 
ergy distribution as that of the star, but follows the spa- 
tial distribution of the dust in the disk. Therefore, the 
observed light curve carries the signature of the lensed 
disk at a small level. Such a light curve is illustrated 
in the top panel of Figure IHI The dashed curve shows 
the classical microlensing light curve of a point source 
(the star) and the solid curve shows light curve of the 



star+disk system, where a fraction of r = 10^^ of the 
star's light is seen in the disk. The light curve shown in 
Figure[Slcan be thought as an extreme version of FigurelHl 
with r ~ 1. The bottom panel of Figure|Blshows the rela- 
tive difference in the magnifications of the star+disk sys- 
tem and the star alone, i.e., the quantity r {Ad/A^, — 1). 
It is essentially at the level of t x 10% near the wings 
and T X 100% near the center of the light curve. The cor- 
responding curve follows similar parameter dependences 
as the ones presented in Figure 0] and can be used in the 
same way to constrain disk properties. 

5. OBSERVATIONAL PROSPECTS 

5.1. Detectability 

In this section we discuss observational requirements 
for the detection of microlensing of circumstellar disks. 
As the ranges of temperature, surface brightness, and 
distance of the systems can be quite large, we will only 
provide some order of magnitude estimates for our fidu- 
cial circumstellar disk located at 8 kpc from us (see Fig^l 
and for a IMq lens located halfway to the Galactic cen- 
ter. Below we will address the detectability of the mi- 
crolensing effects and will consider as an example the cur- 
rent specifications of the upcoming James Webb Space 
Telescope (JWST), which is a 6-meter telescope giving 
access to a wavelength range of 0.6 to 25.5 /im. 

At the Galactic center, a circumstellar disk with radius 
as large as 100 AU has an apparent size of about 0.01", 
which is much smaller than the size of the telescope point 
spread function (PSF). As we have discussed, the prop- 
erties of the disk can be probed with microlensing in 
several ways: from its thermal emission (peaking around 
50/zm) and from the reflected light of the star (peaking 
around 1/im). We will investigate these two possibilities 
separately. 

The JWST cameras allow to observe in the 0.6-25.5/ini 
range. Between 10 and 25/im, the main noise contri- 
butions for observing distant circumstellar disks orig- 
inate from the detector self-emission and the zodiacal 
light (both are more important than the noise level of 
the star). As these two contributions scale linearly with 
the exposure time foxp , the detection signal-to- noise ratio 
(S/N) of the disk is therefore proportional to ^yt^■ In 
this wavelength range, it should be noted that the sensi- 
tivity of the detector decreases exponentially as a func- 
tion of wavelength. Because of this dramatic change, sys- 
tems with warmer disks peaking at shorter wavelengths 
will be much easier to observe. 

Using the sensitivity and noise contributions given 
by the JWST Exposure Time Calculator for the Mid 
Infrared Instrument (MIRI) we have estimated several 
quantities: (i) The exposure time required to detect the 
presence of the disk, i.e., to obtain a 3a detection of the 
infrared excess at 15 and 25.5 /im, respectively. This 
is shown in the bottom panel of Figure |3 (gray lines) 
as a function of the disk-to-star luminosity ratio r. (ii) 
The exposure time required to have a 3cr detection of 
a 10% change in the disk flux, i.e., to detect the mi- 
crolensing signal. This is shown in solid black lines. For 
10""^ < r < 10^^, the exposure time ranges between 1 
and too hours. At 10/im and below, the intensity of the 
blackbody spectrum of the disk decreases rapidly and 
does not allow any detection within a reasonable expo- 
sure time. Finally, we also point out that the current 
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Fig. 7. — The top panel shows, as a function of the disk-to- 
star luminosity ratio r, the quantity age/10"'^''yr, corresponding 
to a rough estimate of the occurrence of microlensing events of 
circumstellar disk systems. The bottom panel presents the required 
exposure time to detect the disk (gray lines) and to detect a 10% 
change in the disk flux (solid black lines) at the Sa level with 
the JWST. The thin and thick solid lines are computed for re- 
emitted light at 15 and 25.5 ^m respectively and include noise 
contributions from the instrument self-emission and the zodiacal 
light. The dashed line corresponds to the microlensing signal seen 
from the reflected light of the disk at 2 fim, for which the main 
source of noise comes from the star light. 



observational limitations come from instrumental noise. 
One of the main noise contributions of MIRI is the in- 
strument's thermal emission. Contrary to the sensitivity 
in optical bands, that of mid/far- infrared detectors is 
expected to keep improving over the years as a result of 
the development of new technologies. Next generation 
mid-infrared detectors might therefore need significantly 
less exposure time than the ones quoted above and may 
even allow us to observe at longer wavelength where the 
emission from disk dominates over that of the star. 

The second way for using microlensing to extract in- 
formation about the disk is to measure the light reflected 
by the disk. In such a case, the signal is the slight devia- 
tion of the microlensing light curve from that of a point 
source, which is typically at the rxlO% level. We as- 
sume that the observation is performed in the K band 
(2.2/im) to minimize the effect of the Galactic extinc- 
tion. At this wavelength, the width of the JWST PSF is 
about 0.1". A ro ugh estim ate based on the Galactic stel- 
lar disk profile in lZhensTet al.l l)20011[ ) shows that, toward 
the Galactic center, the probability that another disk star 
falls within the PSF is less than 0.2%. Inclusion of buldge 
stars would increase the probability by rougly an order 
of magnitude. Therefore, in the case of the lensing of the 
circumstellar disk in reflected light observed by JWST- 
like space telescopes, the main contribution to the noise 
originates from the light of the central star (neglecting 
instrumental and calibration uncertainties). The dashed 
line in Figure [3 shows the exposure time required to have 
a 3(7 detection of a 10% change (caused by microlensing) 
in the flux of the reflected light. For 10"^ < t < 10"^, 
the exposure time ranges between 2 x 10~^ and 2 hours. 



The duration of the microlensing event for a circum- 
stellar disk is of the order of a year. Over such a time 
scale, a few exposures would allow us to sample the disk 
light curve with enough accuracy to constrain a num- 
ber of parameters of the disk like its inner and outer 
radii. We note that the numbers used in this estimate 
are based on a typical circumstellar disk surface bright- 
ness. Brighter disks occur and would therefore reduce 
the exposure time mentioned above. Moreover we have 
based our estimates on a uniform disk. A steeper surface 
brightness profile will have a similar effect as reducing 
the effective size of the disk and therefore increase the 
amplitude of the magnification. 

5.2. Event Occurrence 

The disk-to-star luminosity ratio r of the circumstel- 
lar disk is expected to decrease as a function of the age 
of the system . Observationa lly, it has been shown that 
T cx age~^'^^ l|Zuckerman fc Song 2004). The top axis of 
Figured labels the age of a system corresponding to the 
value of T on the bottom axis. We see that the bright 
phase of the circumstellar disk typically lasts for a few 
to ~100 Myr, therefore the fraction of stars that have 
bright circumstellar disks is determined by the fraction 
of young stars, which is closely related to the star forma- 
tion history. 

For a constant star formation rate over 10 Gyr, the 
fraction of young stars is given by the ratio of the 
stars' age and 10 Gyr. this fraction can be regarded 
as an order of magnitude estimate of the fraction of mi- 
crolensing events that could be those of star-|-disk sys- 
tem, i.e., the occurrence of microlensing events of cir- 
cumstellar disk. This is shown in the top panel of Fig- 
ure [3 The values are only indicative since the star for- 
mation history of the Milky Way is much more com- 
plex than a pure continuous star formation. In the disk, 
there seems to be 5 episodes of enhanced star forma- 
tion over the last 2 Gyr that could be related to den- 
sity waves_aiid/or_2erigaIa£t^^ 

ies llde La Fuente Marcos fc de La Fuente Marcosl2004j) . 
Most young stars stay near the plane where they were 
born since they do not have enough time to diffuse away. 
Therefore, along lines of sight in the Galactic mid-plane 
away from the bulge, the young star fraction is expected 
to be much higher than that shown in the top panel of 
Figure [7| The inner part of the Galaxy is thought to be 
dominated by old bulge stars. This may be true in the 
outer bulge. How ever, clouds and star forming re- 
gions are present l)Banialll977t iLiszt fc Burtonl[T978|) in 
the Galactic center. In the central few hundred parsecs 
of the Galaxy, there are several lines of evidence showing 
rece nt (e.g., 200 Myr ago ) and ongoing star formation 
(see Ivan Loon et al.ll2003l and references therein). So we 
also expect a larger occurrence rate of systems with cir- 
cumstellar disks for lines of sight toward the very center 
of the Galaxy. Taking these recent star formation his- 
tories into account could increase the plotted occurrence 
in Figure |3 by approximately an order of magnitude. A 
more detailed study is necessary for an accurate determi- 
nation of the fraction of young stars along different lines 
of sight. 

Presently several thousands of microlensing events 
have been detected toward the Galactic bul ge by 
OGLE (e.g., lUdalski et all l200?i IWozniak et 311120011 
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ISumi et alJl20p l). MOA (e. g.. [Sumi et al.ll2(ffll . MA- 

CHO (e.g-nr homas eT^ I2nn4j) . and EROS (e.g ., 
lAfonso et al.l2003hV MOA^, OGLE-III^ llUdalskill200.'^ . 

and future experiment will provide us with a large 
amount of data to select candidate events of microlens- 
ing of circumstellar disks to be followed up in near- and 
mid-infrared. The fields monitored by current microlens- 
ing projects are along the lines of sight of low extinction 
regions (e.g., Baade's window), which are typically a few 
degrees away from the Galactic center. To enhance the 
fraction of microlensing events of young stars, surveys 
toward the very center of the Galaxy are desired because 
both the Galactic center and the mid-plane of the disk 
are expected to have higher fraction of young stars. To 
minimize the extinction, such su rveys needs to be done 
in the K band, as advocated bv IGouldl l)1995|) . In ad- 
dition to the fraction of young stars, the microlensing 
optical depths is also enhanced in these surveys. For 
example, the microlensing optical depth of the Galactic 
disk is about 3 times larger toward the Galactic center 
than toward Baade's window (, Gould 1995,) ■ With the 
improvements of infr ared detectors (such a s mosaics of 
infrared arrays, e.g. , iFinger fc Beletid 120031: iDorn et all 
l2004t iMcLeanI 12004(1 . ground-based infrared monitoring 
programs toward the very center of the Galaxy will be- 
come possible in the near future. Space-based infrared 
monitoring programs might also start in a near future 
wit h for example t he Microlensing Planet Finder (MPF, 
see lBennett et aLl l2004'l. They will greatly improve the 
search for microlensing of circumstellar disks. 

Given the different timescales involved in the star and 
disk lensing, it is possible to observe the magnification of 
the disk several months or even a year after the optical 
microlensing event of the star. In order to optimize the 
use of telescope time, the following observational strat- 
egy should be followed: (i) detecting a microlensing event 
in the optical or infrared; (ii) if the star is likely to reside 
in a star-forming region or an II2 cloud that indicates a 
young age, then following up in the near-infrared (mid- 
infrared) in order to detect the presences of a circum- 
stellar disk in reflected (re-emitted) light; (iii) if a disk 
is present, then performing longer exposures to sparsely 
(e.g., once a month) sample the light curve and obtaining 
constraints on the disk properties. In these three steps, 
the second one is crucial. We want the rate of successful 
selection to be as high as possible. For this purpose, we 
can check the (projected) location of the source star to 
see whether it coincides with that of an H2 clouds. We 
can also check the infrared color-magnitude diagram of 
neighbor stars of the source star to see whether there 
is any hint of young stars (e.g., the existence of evolved 
massive stars), (iv) In addition, one can increase the 
probability of observing a microlensing event of a disk by 
preselecting events caused by massive and nearby lenses 
(e.g., by choosing long duration microlensing events) in 
order to have larger projected Einstein ring radii and 
hence higher disk magnifications. 

6. SUMMARY AND DISCUSSION 

We have investigated the microlensing effects on a 
source star surrounded by a circumstellar disk, as a func- 

^ http : / / www, massey. ac . nz / ~ iabond / alert / alert .html 
^ |littp://www.astrouw.edu.pl/~ogle/ogle3/ews/ews.litml| 
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tion of wavelength: 

• In the mid- and far-infrared, where the emission 
of the system is dominated by the cold dusty disk, 
the microlensing light curve encodes the geometry 
and surface brightness profile of the disk. For a 
system located at the Galactic center, typical disk 
magnifications reach 1.1 — 1.2 and the events 
last over ~ 1 year, i.e., much longer than that of 
the star-only case. 

• At around 20/im, where the emission for the star 
and the disk are in general comparable, the differ- 
ence in the emission areas results in a chromatic 
microlensing event. 

• Finally, in the near-infrared and visible, where the 
emission of the star dominates, the small fraction 
of star light directly reflected by the disk causes the 
light curve of the system to slightly deviate from 
that of a point source. In this case, typical de- 
partures from the classical microlensing light curve 
of a point source are at a level of t x 10%, where 
T is the disk-to-star luminosity ratio. For young 
disks, a photometric accuracy of 10"'^ will allow us 
to probe the disk properties. 

We have shown how the disk light curves depend on a 
number of disk parameters like the size, the orientation, 
the inclination, the radius of the inner hole, etc., and have 
pointed out that in some configurations the magnification 
can be substantially higher. The most favorable case is 
the lensing of a small, highly inclined disk with a steep 
surface brightness profile. 

Among the detected microlensing events, there is a 
significant fraction of events with binary lens (e.g., 
■Alcock et al. 2000: Jaros zyhski 2002; Jaroszviisk i et alJ 
l2004|) . In the case of a binary lens, if the separation of 
the binary is much smaller than the circumstellar disk, 
the light curve would largely remain the same as that of 
a single lens. In addition it will better resolve the fine 
features of the disk such as gaps and hot spots, if they 
happen to come close to the caustics of the binary lens. 

The duality between the infrared and optical regimes 
and the long timescale relevant for the disk microlens- 
ing allow us to maximize the observation efficiency by 
only following up systems that have already experienced 
a substantial magnification in the optical. The fraction 
of stars exhibiting a bright circumstellar disk that can be 
detected through the lensing effect is expected to be of 
order a few times 10"'^ with large variations among differ- 
ent lines of sight. A K band micro- lensing survey toward 
the very center of the Galaxy is desired for enhancements 
in both microlensing optical depth and occurrence rate 
of lensing events of circumstellar disks. 

We have investigated the detectability of such events 
with the upcoming JWST telescope. For young disks, 
a few hours of exposure every month would allow us to 
sample the microlensing light curve of the disk and ex- 
tract some disk parameters. In addition, the relatively 
deep exposure of the fields resulting from these observa- 
tions can be used, as a byproduct, for various studies of 
objects toward the Galactic center. 

Microlensing enables us to extend the observation of 
circumstellar disks to a distance far beyond the solar 
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neighborhood, in regions of different environments and 
metalhcities. It has the potential to reveal the diversity 
of circumstellar disks. Studying the amount and nature 
of their dust and their dependence on the environment 
would improve our understanding of the formation and 
evolution of the circumstellar disks and planets. 

We thank Scott Gaudi, Andy Gould, and Bohdan 
Paczyhski for stimulating discussions and useful sugges- 
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